Stadium Billiard with Moving Walls 
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We study the evolution of the energy distribution for a stadium with moving walls. We consider 
one period driving cycle, which is characterized by an amplitude A and wall velocity V . This 
evolving energy distribution has both "parametric" and "stochastic" components. The latter are 
important for the theory of quantum irreversibility and dissipation in driven mesoscopic devices 
(eg in the context of quantum computation). For extremely slow wall velocity V the spreading 
mechanism is dominated by transitions between neighboring levels, while for larger (non-adiabatic) 
velocities the spreading mechanism has both perturbative and non-perturbative features. We 
present, for the first time, a numerical study which is aimed in identifying the latter features. A 
procedure is developed for the determination of the various V regimes. The possible implications 
on linear response theory are discussed. 



I. INTRODUCTION 

Consider the problem of particle in a box, where some 
piece of the wall is deformed periodically in time. As an 
example one may think of a particle in a cylinder with 
a moving piston. The particle has mass m and kinetic 
energy E. The piston is pushed back and forth. The 
velocity in which the piston is displaced is ±1/, and the 
maximum displacement is A. At the end of each cycle 
the piston is back in its original location. 

In the present Paper we are going to study what hap- 
pens to a quantum mechanical particle in such a box, 
during one cycle of the driving. We assume that the par- 
ticle is initially prepared in an energy eigenstate of the 
(un-deformed) box. We shall explain that it is important 
to specify whether the maximum displacement A is less 
or larger compared with De-Broglie wavelength. 

The problem of particlc-in-a-box with a moving wall 
is a prototype example for study of driven systems 
that are described by a Hamiltonian TL{x{t)), where x{ 
is a time dependent parameter. In case of a piston, the 
parameter x is the position of the piston, and we use the 
notation V = \x\. 

There is a lot of interest today in studies of quantum 
irreversibility. This is a relevant issue to the design of 
quantum computers, where the fidelity of a driving cycle 
is important. If at the end of a driving cycle the system 
is back in its initial state, then we say that the driving 
cycle has high fidelity. Obviously, the piston model can 
serve as a prototype example for studying of quantum 
irreversibility [Q. 

The study of one-cycle driving also constitutes a bridge 
to the study of the response to multi-cycle (periodic) 
driving. It should be clear that the long time behavior of 
driven systems is determined by the short time dynamics. 
Therefore it is essential to have a good understanding of 
the latter. 

At first sight one may think that it is simplest to study 
the one-dimensional (ID) box case, also known as 'infi- 
nite square well potential with moving wall' 0, H. The 



case of periodic driving is also known as the 'Fermi accel- 
eration problem' In a second sight one find out that 
the ID-box case is actually the most complicated one. 
As in the case of the kicked rotator (standard map) 
there is a complicated route to chaos and stochasticity. 

Things become much simpler if the motion in the box 
is chaotic to begin with. Driven chaotic systems exhibit, 
as a result of the driving, stochastic energy spreading of 
relatively simple nature . This is the case that we want 
to consider in this Paper. Hence we consider the simplest 
"chaotic box" , which is a two dimensional (2D) "billiard" 
model. Specifically, we are going to introduce a detailed 
numerical study of the one-pulse response of a stadium 
billiard to a shape deformation. The stadium billiard is a 
recognized prototype system for "quantum chaos" stud- 
ies. Our simulations are feasible thanks to a new pow- 
erful technique for finding clusters of billiard eigenstates 
p3| , [T^ . Previous appHcations of this technique, to 
the study of restricted aspects of the present problem, 
have been reported in and in pTf . 

The work which is presented in this Paper, is the first 
numerical study which is aimed in presenting a system- 
atic analysis of non-perturbative features of a time de- 
pendent spreading process |10, |l2|. A procedure 



is developed for the determination of the various time 
stages in the evolution of the energy distribution. This 
allows the identification of the various V regime ("adi- 
abatic", "perturbative", "non-perturbative", and "sud- 
den" ) which were predicted in past theoretical studies. 



II. OUTLINE 

In Sections 3 and Section 4 we define the model 
system, and briefly describe the classical picture. In par- 
ticular we explain that the time dependent features that 
we study in this Paper are purely quantum-mechanical, 
and not of semiclassical origin. 

In Section 5 we define the main object of our study, 
which is the energy spreading kernel Pt (n|m) . See Eq. (0) . 
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Regarded as a function of the level index n, it is the en- 
ergy distribution after time t, while m is the initial level. 
We also define the sqrt of the variance 5E{t), and the 
50% probability width r{t), that characterize the evolv- 
ing energy distribution. 

Our three phase strategy for analysis of energy spread- 
ing is presented in Section 6. The three phases are: 

(I) Study of the band profile. 
(II) Study of parametric evolution. 
(Ill) Study of the actual time evolution. 

The relevant information regarding " Phase I" is summa- 
rized in Section 7. Various approximations for Pt{n\'m), 
and in particular the notion of "parametric evolution" 
are presented in Section 8. The relevant information 
regarding "Phase 11" is summarized in Section 9. The 
main concern of this Paper is with "Phase III". 

For a given V one should be able to characterize 
the nature of the dynamical scenario. The theoreti- 
cal considerations regarding this issue have been dis- 
cussed in m ^ |ll], |l^ , and for a concise review see . 
Rather than duplicating these discussions, we are going 
to present in Section 10 a phenomcnological definition, 
and a practical procedure, for the identification of the 
V regimes. We would like to emphasize that this is the 
first time that the different dynamical scenarios (corre- 
sponding to the different V regimes) are illustrated in a 
numerical simulation. The only other numerical studies 
(mainly |l^) were too restricted in scope, and did 
not contain analysis of the stages in the evolution of the 
energy distribution. 

Section 11 and Section 12 question the applicability 
of linear response theory (LRT) to the analysis of the en- 
ergy spreading. Further discussion of "non perturbative 
response" , and conclusions, are presented in Section 13. 



III. THE SEMICLASSICAL PICTURE 

Consider a classical particle inside a box. Its kinetic 
energy is E and the corresponding velocity is ve = 
\/2Elm. The shape of the box is externally controlled. 
The control parameter is denoted by x. We assume that 
X has units of length, such that V — \x\ is the typical wall 
velocity. Obviously different parts of the wall may have 
different velocities. In case of the "piston model" only 
one piece of wall is moving (either inward or outward). 
However, we are not interested in the trivial conservative 
work which is being done, but only in the irreversible 
work. Therefore, rather than analyzing an actual "pis- 
ton model" configuration, it is wiser to consider a volume 
preserving deformation. In such case the conservative 
work is zero, and the major issue is the irreversible effect 
which is explained below. 

As a results of the collisions of the particle with the 
deforming walls, there is a stochastic-like diffusion in en- 
ergy space. The explanation is as follows: Each time 



that the particle collides with the moving wall it either 
gains or looses energy. To simplify the presentation let 
us assume head-on collisions. The change in energy is 
±2mvEV depending on whether the wall is moving in- 
ward or outward at the point of the collision. For volume 
preserving deformation the ergodic average over SE gives 
zero. Thus we have random- walk in energy space where 
the steps are ±2mvEV . This leads to diffusion in energy 
space. As explained in this diffusion is biased (the 
diffusion is stronger for larger E). This leads, in the long 
run, to a systematic (irreversible) increase of the average 
energy. For a more detailed presentation that takes the 
box geometry into account see Q . 

Do we have a corresponding (semiclassical) picture in 
the quantum-mechanical case? Again we remind the 
reader that we assume a volume preserving deformation. 
This implies that the energy levels of the system do not 
have a collective "upward" or downward" change as a 
result of the deformation. The physics in which we are 
interested is related to the transitions between different 
levels. The question that we ask is whether these transi- 
tions are " classical-like" . 

Let us assume that we start with an eigenstate whose 
energy is E. Semiclassically it is as if we start with a 
microcanonical preparation. If the dynamics is of classi- 
cal nature, then we expect that after a short time some 
of the probability will make a transition E ^ E' such 
that \E' — E\ ^ 2mvEV. Naturally such description is 
meaningful only if the energy scale 2mvEV is much larger 
than the mean level spacing. This leads for ID-box to the 
condition 

V > h/mL (1) 

In the general case (eg 2D box), having 2mvEV much 
larger than the mean level spacing is not a sufficient con- 
dition for getting semiclassical behavior. Still, non-trivial 
analysis |22| reveals that the condition for getting semi- 
classical spreading in the general case is the same as in 
the ID case, namely given by Eq.(|l|). 

The above Eq.(|^) is a necessary condition for hav- 
ing semiclassical transitions. In order to actually vuit- 
ness semiclassical transitions it is also necessary to have 
a time period much larger compared with the ballistic 
time {A/V ^ L/ve), and to have an amplitude that 
is much larger compared with De-Broglie wavelength 
{A ^ h/{mvE))- These additional conditions can be 
satisfied only in the semiclassical regime which is defined 
by Eq.(^, else they are not compatible. 

IV. THE NUMERICAL MODEL 

As a specific example for chaotic box, we consider the 
quarter stadium billiard. We define x as the length of the 
straight edge, and adjust the radius parameter such that 
the total area is kept constant. For numerical reasons 
we do not analyze this model "literally" but rather con- 
sider, as in previous study pq|, a linearized version of the 
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quarter stadium billiard Hamiltonian. Namely, we study 
a model Hamiltonian that has the matrix representation 



practical purpose, given f(t), it is convenient to associate 
with it the a spectral function 



n{x{t)) ^'E + Sx{t)F 



(2) 



Here E is an ordered diagonal matrix, that consists of 
the eigen-energies En of the quarter stadium billiard with 
straight edge x = 1. The eigen-energies were determined 
numerically. The perturbation due to Sx deformation, is 
represented by the matrix F. Also this matrix has been 
determined numerically as explained in [ p^ . We note 
that the fingerprints of the classical chaos are present in 
the statistical properties of the matrices E (level statis- 
tics) and F (band structure). The latter is discussed in 
Section 7. 

The linearization Eq.(^ of the billiard Hamiltonian 
can be regarded as a valid approximation if the wall 
displacement parameter 6x is small compared with De- 
Broglie wavelength. This automatically excludes the pos- 
sibility of addressing the semiclassical regime which has 
been discussed in the previous section. In our simulation 
the De-Broglie wavelength is roughly 0.1. This implies 
that (at best) the maximum driving amplitude that can 
be allowed is ^ = 0.2, so as to have |i5a;(i)l < O-l- [Note 
that for presentation purpose we later re-define 5x{t) as 
Sx{t) — Sx{Q).] In the simulations we indeed have A = 0.2, 
meaning that we allow a deformation which is compara- 
ble or possibly somewhat larger than that "allowed" by 
the linearization. 

However, inspite of the fact that Eq.(||) may be a 
"bad" or even "inadequate" approximation to real- world 
physics, it is still a totally "legitimate" Hamiltonian from 
mathematical point of view. Moreover, the numerical 
model Eq.(^ contains all the physical ingredients that 
are relevant for the aim of the present study. 

To summarize: In this paper we consider, as far as for- 
mulation is concerned, a chaotic billiard driven by volume 
preserving shape deformation. On the other hand, as far 
as numerics is concerned, we analyze a specific quantum 
mechanical model which is defined by Eq.(^). The nu- 
merical model is motivated by the quarter stadium bil- 
liard system Hamiltonian, but still it is not literally the 
same model. In particular, Eq.(H) does not possess the 
the semiclassical regime which has been discussed in the 
previous section. 



V. THE EVOLVING ENERGY DISTRIBUTION 

In case of time independent Hamiltonian, the energy 
distribution does not change with time. In order to have 
an "evolving" energy distribution, we have to make Sx{t) 
time dependent. One possibility is to assume linear driv- 
ing. In such case we write Sx{t) — Vt. But more gen- 
erally, for a cycle, we write Sx{t) ~ A x f{t), with the 
convention /(O) = f{T) = 0. In some equations below, 
whenever a linear driving is concerned, A can be replaced 
by V, which assumes the particular choice f(t) = t. For 



f{t')e"^* dt' 



(3) 



Some useful cases are summarized in Appendix A. In this 
Paper we are p rimarily interested in the case of triangular 
pulse Eq.dAlOD. 

Given the model Hamiltonian Eq.(|^), and the driving 
scheme 5x{t), we can calculate the unitary evolution op- 
erator U{t). This kernel propagates " wavefunctions" in 
time. Equivalently, we can describe the quantum me- 
chanical state using a probability density matrix. The 
propagator of the latter is denoted by A(f). In order 
to describe the evolution of the energy distribution we 
define the kernel 



Pt{n\m) = trace(p„A(i)pm) 

= |(n(x(t))|C/(t)|m(x(0)))p 



(4) 



where p„ can be interpreted as either the probability ma- 
trix or as the corresponding Wigner function that repre- 
sents the eigenstate \n). In the latter case the trace oper- 
ation should be interpreted as dPdQ / {2'KhY integration 
over phase space, where (Q, P) are the canonical coordi- 
nates of the particle, and d = 2 is the dimensionality. 

We shall use the notation Pt{r) = Pt{n — m) = 
Pt{n\m), with implicit average over the reference state 
TO. We shall refer to Pt(r) as the "average spreading 
profile" . Whenever we have a wide distribution, we dis- 
regard the distinction between "energy difference" and 
" level difference" , and make the identification 



En — Eni = huJn 



A X 



(5) 
(In 



where A is the mean level spacing, and r = n — 
our numerics A w 7.22). 

Fig.l displays the evolution as a function of time. 
Each column is a profile of the probability distribution 
Pt(r). The first 7 time steps are log spaced, while the rest 
are linearly spaced. The V = 100 evolution is predom- 
inantly parametric. As V becomes smaller and smaller 
the deviation from parametric evolution becomes larger 
and larger. Some representative spreading profiles are 
presented in Fig. 2. 

There are various practical possibilities available for 
the characterization of the distribution Pt{r). It turns 
out that the major features of this distribution are cap- 
tured by the following three measures: 

V{t) = Pt(r = 0) (6) 
r{t) = 50% probabihty width (7) 

1/2 

(8) 



SE{t) = Ax ^r^Ptir)^ 



The first measure is the survival probability V{t). It 
is mentioned here just for completeness of presentation. 
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The second measure r(t) is the energy width of the cen- 
tral r region that contains 50% of the probabihty. It is 
calculated as T{t) = (r75%-r25%)A, where r25% and ri5% 
are the values for which the cumulative energy distribu- 
tion equals 25% and 75% repectively. Finally the energy 
spreading SE{t) is defined above as the square-root of the 
variance. 

Fig. 3 shows how the width T{t) of the profile evolves. 
The lower panel of Fig. 3 is a log-log plot. We see clearly 
that up to 6xc = 0.006 the width is one level, which 
is an indication for the applicability of standard first- 
order perturbation theory. For larger Sx several levels 
are mixed non-perturbatively, and therefore the width 
becomes larger than one. 

Fig. 4 shows how the spreading SE{t) of the profile 
evolves. The lower panel displays the relative spreading, 
which is defined as the ratio between the actual spreading 
and the parametric one. The notion of "parametric evo- 
lution" is defined in the next paragraph. As V becomes 
smaller, the departure from the parametric behavior hap- 
pens earlier. 

The "sudden limit" {V — > oo) of Pt{n\m) will be de- 
noted by P{n\m). We shall refer to it as the "paramet- 
ric kernel" . It is formally obtained by the replacement 
A{t) I— > 1, or equivalently J7(t) 1 in Eq.(^, namely, 

P{n\m) = trace(p„pm) 

= |(n(:r(t))|m(x(0)))p (9) 

Obviously its dependence on t is exclusively via Sx — 
x{t)—x{0), irrespective of V. The parametric evolution of 
P(r) versus Sx for a deforming billiard has been studied 
in . The numerics in this Paper can be regarded as the 
extension of the numerics of [0 to the case of finite V. 

VI. THREE PHASE STRATEGY 

In the future we want to have a theory that allows the 
prediction of the (numerically simulated) time evolution. 
The minimum input which is required for such theory is 
the band profile of the F matrix. For precise definition of 
the bandprofile see Appendix B. Note the existence of a 
very efficient semiclassical recipe to find the band profile. 
This can save the need for a tedious quantum mechanical 
calculations. 

Using the bandprofile one hopes to be able, in the 
"second phase" to calculated the parametric kernel 
P{n\m). The bandprofile does not contain information 
about the correlations between the off-diagonal elements. 
Therefore one has to make the so-called RMT conjecture, 
namely to assume that the off-diagonal elements are ef- 
fectively uncorrelated. It turns out that, upon using such 
conjecture, the non-universal features of the parametric 
evolution are lost. Still, the obtained results are quali- 
tatively correct, and therefore such an approach is legit- 
imate as an approximation. 

However, finding the parametric kernel P(n\m), using 
the bandprofile as an input, is not the subject of this 



paper. Therefore, as a matter of strategy, we would like 
to take the numerically determined parametric evolution 
as an input. Given the parametric P(n\m), the question 
is whether we can calculate the actual evolution Pt{n\m) 
for any finite V. 

In previous works (see review in ) we gave a negative 
answer to the above question. We have claimed that non- 
perturbative features of the dynamics cannot be deduced 
from the parametric analysis. To explain this observation 
let us assume that we have two model Hamiltonians, say 
'Wphyi,icai(a;) and 7iartificiai(a^)- Let us assume further that 
the two models have the same bandprofile and the same 
parametric kernel P(n\m), Still we claim that for finite 
V the two Hamiltonians may generate different temporal 
kernels Pt{n\m). In particular, it has been argued that 
this is in fact the case if Tiartiticiai (a;) is an effective RMT 
model which is associated with 7iphyBicai(a;)- 

In past publications the manifestation of non- 
perturbative features in case of driven physical (non- 
RMT) models has not been investigated numerically. On 
the theoretical side it is an open subject for further re- 
search The purpose of the following sections is 
to present a strategy for analysis of numerical simula- 
tions, that paves the way towards a theory for the non- 
perturbative aspects of the energy spreading. 



VII. THE BANDPROFILE 

The bandprofile of F is described by the spectral func- 
tion C{lu), which is a Fourier transform of a correlation 
function C(t). See Appendix B. If the collisions are 
uncorrelated, as in the case where the deformation in- 
volves only a small surface element, then C(r) is a delta 
function, and the band profile is fiat. This is not the 
case in the present model. The correlations between col- 
lisions cannot be neglected, and therefore C(r) equals 
delta function plus a smooth component Due to 

the smooth component the band profile has a very pro- 
nounced non-universal structure. See Fig. 5. 

For large enough co the contribution of the non- 
universal component vanishes. Thus the bandprofile 
should possess fiat tails. These tails reflect the presence 
of the delta- function component in C(t). For numerical 
reasons, due to truncation, the bandprofile of our sys- 
tem is multiplied by a Gaussian envelope. Therefore the 
tails of the effective bandprofile are not flat, but rather 
vanishingly small. The lack of fiat tails in the numerical 
model can be loosely interpreted as having 'soft' rather 
than 'hard' walls. 

Thus our system is characterized by a finite bandwidth. 
This is actually the generic case. Namely, for generic 
("smooth") Hamiltonian the correlation function C(r) is 
non-singular, which implies finite bandwidth. 
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VIII. APPROXIMATIONS 

A major issue in the studies of energy spreading is 
the knowledge how to combine tools or approximations 
in order to understand or calculate the kernel Pt(n|rn). 
In particular we have the following approximations: 

• The perturbative kernel P^'^{n\m). 

• The semiclassical kernel Pp\n\m). 

• The Gaussian kernel P^*°{n\m). 

• The parametric kernel P(n\m). 

The parametric kernel P{n\ni), that corresponds to the 
sudden limit {V oo) has already been defined in 
Eq.(^. The other kernels will be defined in the present 
section. Obviously the kernels listed above are very dif- 
ferent. Our aim is to clarify in what regime (V), in what 
time stage (t), and in what energy region (r), which of 
them is the valid approximation. 

The semiclassical kernel P('''(n|m) is defined and ob- 
tained by assuming in Eq.(|) that the Wi gner functions 
can be approximated by smeared microcanonical distri- 
butions. See Ref.^, |l^, |3| for further details and nu- 
merical examples. Whenever Pt{n\m) ^ P^'''-{n\m) we 
say that the spreading profile is " semiclassical" . In order 
to have a valid semiclassical approximation in case of bil- 
liard systems the displacement Sx of the walls should be 
much larger than De-Broglie wavelength [0. Thus, for 
reasons that were explained in Section 3, the semiclas- 
sical approximation is not applicable for the numerical 
model that we consider in this Paper. 

The perturbative kernel P['''{n\m) is obtained from 
perturbation theory The defining expression is: 



Prin\m) 



t \ 



(10) 



where T is determined by normalization. If we make the 
replacement F i-^- we get the standard result of first 
order perturbation theory (see Eq.(^5|)). The presence 
of r in the denominator reflects corrections to infinite 
order. We note that P/"'(n|TO) constitutes a generaliza- 
tion of Wigner Lorentzian. We indeed would get from 
it a Lorentzian if the bandprofile were flat with no finite 
bandwidth. 

Whenever Pt{n\m) ~ P/"'(n|TO) we say that the 
spreading profile is "perturbative". The perturbative 
structure is characterized by having separation of scales: 



r(t) < SE{t) 



(11) 



We say that Pt{n\m) has a "standard" perturbative 
structure Eq. ( p^ ) , which is given by first order perturba- 
tion theory, if F < A. This means that more than 50% 
probability is concentrated in the initial level. We use 
the term "core-tail" structure if we want to emphasize 
the existence of a finite non-perturbative "core" region 



|r| < F/A. The core width F, as determined by "pertur- 
bation theory" by imposing normalization on Eq.(|o|), 
constitutes a rough estimate for the width F(t) of the 
energy distribution (Eq.(^). 

The spreading SE{t) for the core-tail structure is deter- 
mined by the tail region, which is defined as \r\ ^ F/A. 
If Eq.(lO) were a Lorentzian, it would imply 5E[t) — oo. 
This is of course not the case, because the spectral func- 
tions provide a physical cutoff. Thus, the core-tail struc- 



ture which is described by Eq.(lO) is characterized by a 
" tail" component that contains a vanishingly small prob- 
ability but still dominates the variance. 

If we do not have the separation of energy scales 
Eq. ( |lTl) , then perturbation theory becomes useless (|, ^ . 
In such case Pt(n|m) becomes purely non-perturbative. 
In order to determine Pt{n\'m) we have to use tools that 
go beyond perturbation theory. In particular, for long 
times one can justify 1^, |l^ a stochastic approximation, 
leading to the Gaussian kernel 



Pr{n\m) 



A 



exp 



1 / E n — Em 

2 \ SE{t) 



(12) 



Note that for very long times this Gaussian should be 
replaced by an appropriate solution of a diffusion equa- 
tion. But this is not relevant to our simulations. For 
a critical discussion that incorporates estimates for the 
relevant time scales see |l^ . 

A final remark: For a Gaussian profile the 50% width 
is T{t) = 1.35 X 5E(t). However, whenever a non- 
perturbative structure is concerned it is better, in order 
to avoid confusion, not to use the notation T(t). The 
notation F has been adopted in the common diagram- 
matic formulation of perturbation theory. In case of non- 
perturbative structure this formulation becomes useless, 
and therefore the significance of F, as a distinct energy 
scale, is lost. 



IX. ANALYSIS OF THE PARAMETRIC 
EVOLUTION, 5x REGIMES 

The parametric evolution of P(n|m) is illustrated in 
Fig. 1 (upper panel) . For very small 5 x we observe clearly 
a standard perturbative profile (Eq.(p^)), whose struc- 
ture is just a reflection of the bandprofile. By defini- 
tion this holds in the standard perturbative (paramet- 
ric) regime 5x < Sxc- From the numerics (looking on 
the lower panel of either Fig. 3 or Fig. 6) we find that 
Sxc = 0.006. 

For Sx > 6xc the initial level starts to mix with neigh- 
boring levels. As a result a non-perturbative 'core' com- 
ponent starts to develop. Thus we obtain a core-tail 
structure. The tail is the perturbative component. The 
main component of the tail is the "first order tail" which 
can be calculated using first order perturbation theory. 
The tails in the vicinity of the core are growing slower, 
which can be regarded as a suppression of core-to-tail 
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transitions due to the mixing. There are also higher or- 
der tails that can be neglected. 

Having 5E ^ F is an indication that 6E is domi- 
nated by the (perturbative) tail component of the core- 
tail structure Eq.(|l0|). The condition T ^ 6E, can be 
re- written as Sx <C Sx^^t, which constitutes a definition 
of the extended perturbative (parametric) regime. Un- 
fortunately, the strong inequality T <C 6E is nowhere 
satisfied in our numerics. Therefore the numerical def- 
inition of Jxprt becomes ambiguous. One may naively 
define Sx^^t by transforming the weak inequality T < SE 
into a weak inequality 6x < Sx^^t. But this procedure is 
numerically meaningless: The quantities T and SE are 
energy scales. As such their definition is arbitrary up to 
a prefactor of order unity. After some thinking one real- 
izes that the only practical definition for Sx^^^t is as the 
6x where the P/"' based prediction of 6E becomes sig- 
nificantly less than the actual spreading. From the nu- 
merics (Fig. 6) we find that dx^^t = 0.05. This definition 
is not ambiguous numerically because we compare "vari- 
ance based measure" (SE) to "variance based measure" 
((5i?''"), rather than comparing "variance based measure" 
(SE) to a "width measure" (F). 

For 6x ^ SXp^t the P^'* based calculation of 6E gives 
saturation. This (non-physical) saturation is the conse- 
quence of having finite bandwidth. It should be regarded 
as an artifact that reflects the limited validity of pertur- 
bation theory. 

In the non-perturbative (parametric) regime, namely 
for 5x > (5a;prt , the tails are no longer the dominant com- 
ponent. This is associated with a structural change in 
the spreading profile. Looking at the upper panel of Fig.l 
we observe a core-tail structure up to the 13th time step. 
Then, the secondary lobe of the bandprofile is swallowed 
by the core. This happening is refiected in Fig. 6 (upper 
panel) by a boost in the width F. [The notion of "sec- 
ondary lobe" should be clear by looking at the \r\ ^ 20 
region in the upper panel of Fig. 2. Unfortunately its vis- 
ibility in the corresponding image (upper panel of Fig.l) 
is quite poor. Still it is possible to follow the evolution of 
the r ~ 20 tail region, and to realize that it is swallowed 
by the core] . 

To have a weak inequality dx > Sx^^^ is not quite the 
same as to be in the (deep) non perturbative regime 
where Sx ^ Sx^^t. For Sx > Sxp^t we do not get a purely 
non-perturbative structure. Inspite of the failure of the 
core-tail picture Eq. ( p^ ) , we still can make a phenomena- 
logical distinction, between" core" and "tail" regions. One 
clearly observes that the tail region is "pushed" outside 
because of the expanding core. This non-perturbative ef- 
fect is not captured by Eq. (|l0|) . In a previous study |l^ 
the crossover from a core-tail structure to a purely non- 
perturbative (semiclassical) structure was quite abrupt. 
In the present study the "deep" non-perturbative regime, 
where all the tail components disappear, is not accessible 
due to numerical limitations. A "purely non perturba- 
tive" structure would be obtained if all the tail compo- 
nents were swallowed by the expanding core. 



X. ANALYSIS OF ACTUAL EVOLUTION, 
V REGIMES 



We start this section with a qualitative description of 
the actual evolution, which is based on looking in illus- 
trations such as in Fig.l and Fig. 2. Later we present the 
quantitative analysis. In general one observes that for 
finite V the evolution acquires stochastic-like features. 
At intermediate times (0 < i < T) the spreading profile 
contains both parametric and stochastic " components" . 
The parametric component is "reversible". In contrast 
to that the stochastic component of the spreading is not 
affected by the velocity reversal, and may lead to a final 
Gaussian distribution. 

Let us describe what happens to the evolution, as a 
function of 6x = Vt, as we make simulations with smaller 
and smaller V. On the basis of Eq.(|l^) we expect to ob- 
serve a modulation of the tails of Pt{r) in a way that is 
implied by the presence of envelope Ft{oj). This modu- 
lation is characterized by the energy scale Ti/t = hV/Sx. 
If we make V further and further smaller, the secondary 
lobes of the tails [see previous section for definition] do 
not have a chance to become visible, because they are 
suppressed by the narrower envelope Ft{uj). For small 
enough V only the main core component of the bandpro- 
file is left visible. For very small V the spreading profile 
becomes very close to a Gaussian shape in a very early 
stage of the evolution. 

The spreading at the end of the pulse period is very 
small both is the sudden limit (very large V , correspond- 
ing to multi-period driving with large frequency), and 
also in the adiabatic limit (very small V). Fig. 2 illus- 
trates representative spreading profiles, which are ob- 
served in the end of the pulse period. Al so d isplayed 
are the first-order perturbative profile (Eq.(D5)), and a 
Gaussian with the same width. Note that the first-order 



perturbative profile Eq.(D5), unlike Eq.(|lO[), does not 
have a proper normalization. We can define a 'differ- 
ence measure' in order to quantify the deviation of the 
(actual) lineshape from Gaussian shape. One possible 
definition is 



/E ^"°(^) [iog(^t(0) - iog(pr(0)]' (13) 



The lower panel of Fig. 8 shows that for y < 20 the 
spreading profile at the end of the pulse is close to Gaus- 
sian shape, while for > 20 this profile is predominantly 
of perturbative nature. 

We should address now the issue of V regimes. Namely, 
for a given V we would like to characterize the nature 
of the dynamical scenario. Below we shall introduce 
a practical procedure for the identification of the vari- 
ous regimes. We are going to explain that in the nu- 
merical simulation we observe four different V regimes: 
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• Adiabatic regime (V^ < 3). 

• Perturbative regime (3 < V < 7). 

• Non-pertm-bative regime {7 < V < 20). 

• Sudden regime (20 < V). 

Before we get into details, we would like to make a con- 
nection with the theoretical discussion of regimes in . 
There we have considered a (sinusoidal) periodic driv- 
ing that is characterized by amplitude A and frequency 
n = 2'k/T. For sinusoidal driving the root-mean square 
"velocity" is = A^/\/2. So fixing A and chang- 
ing V in the present paper, is completely analogous to 
fixing A and changing VL in Ref . p| . Thus, the four 
V regimes listed above correspond to a horizontal cut 
A = const > ^prt) in the (il,A) regime diagram of 



Ll| (see Fig. 1 there). One should realize that the "per- 
turbative regime" is in fact the "linear response (Kubo) 
regime" as defined in , and that the " sudden regime" 
corresponds to the regime of vanishing (high frequency) 
response. 

In the adiabatic V regime the spreading is due to near 
neighbor level transitions, and disregarding extremely 
short times (for which we ca n ap ply standard first or- 
der perturbation theory (Eq.(|D5|)), it looks stochastic. 
This means that Eq. (^ is a quite satisfying approxima- 
tion. For further details regarding the identification of 
the adiabatic regime see Section 12. 

Outside of the adiabatic regime there is a time stage 
where perturbation theory (Eq.(p^) is valid. But after 
this time we have to use theoretical considerations that go 
beyond perturbation theory. In the following paragraph 
we divide the non-adiabatic regime into "perturbative" 
V regime, "non-perturbative" V regime, and "sudden" 
regime. The basis for this distinction is the timing of 
the departure from parametric evolution. The departure 
from parametric evolution is best illustrated by the lower 
panel of Fig. 4. 

The perturbative V regime is defined by the require- 
ment of having the departure from parametric evolution 
happen before the breakdown of Eq. (|lO|) . Consequently, 
in this regime the departure time can be deduced from 
Eq.(|lO|). For example, let us assume a simple bandpro- 
file. In such case departure time is just the inverse of the 
bandwidth, which implies via Eq.( p33| ) that it is simply 
the classical correlation time t^i of the chaotic motion. In 
the actual numerical analysis the bandprofile is not " sim- 
ple" , but has some structure. Rather than debating over 
the definition of Td , it is more practical to determine the 
departure time by inspection of Fig. 4. For convenience 
we have indicated the location of 5x-p^t by a vertical line. 
For V < 7 the departure from parametric evolution hap- 
pens before this line. This way one can determine the 
upper border {V — 7) of the perturbative regime. 

The non-perturbative V regime is defined by the re- 
quirement of having the departure from parametric evo- 
lution happen after the breakdown of Eq.(lO). In other 
words, it means that this departure cannot be captured 



by perturbation theory. A more careful definition of the 
non-perturbative regime should take into account the 
driving reversal at t ^ T/2. As explained in the next 
paragraph, the non-perturbative regime is further re- 
stricted by the requirement of having the departure from 
parametric evolution happens before the driving reversal. 

The distinction between "perturbative" and "non- 
perturbative" V regimes is related to the possibility to 
"capture" the stochastic features of the energy spreading 
by perturbation theory. Non-perturbative features of the 
energy distribution during intermediate times are of no 
relevance if they are of parametric (non-stochastic) ori- 
gin. Let us look for example, in the upper panel of Fig.l. 
We clearly see that the spreading profile at the end of 
the pulse (at t = T) is of standard perturbative nature. 
All the non-perturbative features that develop during the 
first half of the driving cycle, are completely reversed in 
the second half of the cycle. 

Conceptually the simplest way to set a criterion for get- 
ting at the end of the pulse a perturbative structure, is 
to use "fixed basis" perturbation theory (see Appendix C 
of ll^). If we adhere to the present (more physical) ap- 
proach of using x-dependent basis, an equivalent method 
Jl2t is to determine the sudden time t^^„. This is the time 
to resolve the expanding " core" . It is determined as the 
time when the following inequality is violated: 



r(<) < h/t 



(14) 



The regime where the condition t^^^ <C T is violated is 
defined as the " sudden regime" . In the sudden regime all 
the non-perturbative features are of parametric nature, 
and therefore the validity of perturbation theory survives 
at the end of the pulse. In the sudden regime the de- 
parture from parametric evolution becomes visible only 
after the driving reversal. One can determine the sud- 
den regime simply by looking at the lower panel of Fig. 8. 
The departure of the energy distribution (at t = T) from 
Gaussian shape is correlated with getting into the sudden 
regime. 

The discussion above of V regimes was based on look- 
ing for the breaktime of perturbation theory Eq.(lO), on 
the one hand, and looking for the departure from para- 
metric behavior on the other hand. The departure from 
parametric behavior is an indication for the appearance 
of a predominant stochastic component in the spread- 
ing profile. This departure does not imply that we 
get a Gaussian line shape. In the perturbative regime 
we get the Gaussian line shape after the breaktime of 
perturbation theory, which happens after the departure 
from parametric behavior. Therefore in the perturba- 
tive regime there are 3 stages in the evolution: a para- 
metric stage, a perturbative stochastic stage, and a gen- 
uine stochastic stage. In contrast to that, in the non- 
perturbative regime we do not have an intermediate " per- 
turbative stochastic stage", because the departure form 
parametric behavior happens after the breakdown of per- 
turbation theory. 
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XI. LRT FORMULA 

The general LRT formula for the variance of the 
spreading is 



A' 



(15) 



The proportionality SE cx A reflects having "linear re- 
sponse" . Two spectral functions are involved: One is 
the the spectral content of the driving (Appendix A), 
and the other is the power spectrum of the fluctuations 
(Appendix B). The latter is the Fourier transform of a 
correlation function C(t). 

A special case is the sudden limit {V = oo) for which 
Ft{ui) = 1 and accordingly 



SE{t) = ^C{t = 0) X A 



(16) 



Another special case is the response for persistent (ei- 
ther linear or periodic) driving where Ft(uj) = tx 2^6(1^) 
implies diffusive behavior: 



SE = ^/2DEt 



(17) 



In such case the expression for I?^; is known as Kubo 
formula, leading to a fluctuation-dissipation relation. 

Finding the conditions for validity of LRT is of major 
importance. This should be regarded as the first step 
in the analysis of the response of a driven system. The 
classical derivation of Eq.(|l5|) is quite simple, and for 
completeness we present it in Appendix C. For sake of 
the following discussion one can assume that the classical 
" slowness" conditions for the validity of classical LRT are 
satisfied. 

The quantum mechanical derivation of this formula is 
much more subtle |, |l0[ |ll|, m and leads to the dis- 
tinction between " adiabatic" and " (extended) perturba- 
tive" and " non-perturbative" regimes. The semiclassical 
limit is contained in the latter regime. It is known that 
the LRT formula does not hold in the adiabatic regime 
[20 1, but it is valid in the (extended) perturbative regime 
[ 10, |l^ . It is not necessarily valid in the non-perturbative 
regime [|l^ |ll[, but if the system has a classical limit, it 
must be valid again in the semiclassical regime. See fur- 
ther discussion in the next sections. 



XII. ANALYSIS: COMPARISON WITH LRT 

Using the bandprofile as an input we can calculate 
the spreading using Eg. (|l5|) . In Fig. 6 the calculation 
is done for the parametric evolution (dashed line given 
by Eq.(|l^)), while in Fig. 7 the calculation is done with 
Eq.(p^ for finite V values. The latter figure should be 
compared with the upper panel of Fig. 4. The agreement 
is very good unless the velocity V is small. See later dis- 
cussion of the QM-adiabatic regime. The spreading 6E 
at the end of the pulse is better illustrated in Fig. 8. 



There are three possible strategies for evaluating the 
bandprofile. The first is to use the semiclassical strategy 
with Eq.(B3). The second is to make a careful numer- 
ical evaluatio n of the matrix elements, and to use the 
definition Eq.(Bl). This gives the thin line in Fig. 5. 



However, the most practical and economical procedure is 
simply to deduce the bandprofile from the spreading pro- 
file P{r) that correspond to the smallest Sx value. This 
is the thick line in Fig. 5, which we regard as the most 
appropriate estimate. 



We found out that the LRT formula Eq.(15) is not sen- 
sitive to the way in which the bandprofile is evaluated, 
unless the |n — m| = 1 elements of the F„m matrix domi- 
nate the result. Let us denote by a the root-mean-square 
magnitude of these matrix elements. The sensitivity to 
a happens in the V regime which is determined by the 
condition V < (A)'^ /{ha), which is the adiabaticity con- 
dition. This sensitivity can be used as a practical tool for 
the determination of the adiabatic regime. In the upper 
panel of Fig. 8 we display the result of the LRT calcula- 
tion upon setting a = 0. We see that the LRT calculation 
implies QM-adiabatic behavior for V < 3. 

If V is well away from the adiabatic regime, then near- 
neighbor level transitions have negligible contribution. In 
such case the LRT calculation is not sensitive to the ex- 
act value of cr. As V becomes smaller, there is a larger 
relative weight to the near-neighbor matrix elements. 

Deep in the QM-adiabatic regime the mean-level en- 
ergy difference is resolved much before the levels are 
mixed, and therefore, as a result of recurrences, the prob- 
ability stays concentrated in the initial level. This is of 
course a leading order description. In fact we cannot 
neglect higher order corrections. 

In the adiabatic regime the spreading is dominated by 
transition between neighboring levels (whereas outside of 
the adiabatic regime the contribution of neighboring level 
transitions can be neglected). Therefore, it is only in the 
adiabatic regime where the quantum-mechanical calcula- 
tion (using Eq.([l^)) gives results that are different from 
the classical expectation [Q. Is it possible to witness a 
regime where the quantum mechanical (rather than the 
classical) LRT prediction is observed? Apparently the 
answer is positive, but not in a typical numerical exper- 
iment. The reason is that almost always the Landau- 
Zener (non LRT) mechanism for energy spreading takes 
over |g. 

We can verify the dominance of Landau- Zener mech- 
anism as follows. The theoretical prediction is 5E(t) cx 
\/DEt with De <x V^l^ . In the upper panel of Fig. 8 
we plot the spreading at the end of half pulse period 
(t = T/2 = A/V) and at the end of full pulse period 
[t ^ T ^ 2A/V). Hence we expect 5E cx V'^''^. The 
agreement with this expectation is quite good if we con- 
sider the spreading after half pulse period. At the end of 
full pulse period we get that the spreading is larger than 
expected. This is apparently due to the non-adiabatic 
nature of the V i— > —V switching. 
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XIII. DISCUSSION AND CONCLUSIONS 



The numerical study that we have presented should be 
regarded as an application of a general procedure for the 
analysis of energy spreading. In the summary below we 
re-order the stages in this analysis in the way which is 
implied by the results of the previous sections. 

The first step in the analysis is to find the band pro- 
file. This can be done using the semiclassical recipe (Ap- 
pendix B) without any need to make heavy numerical 
simulations. Then it is possible to determine SE{t) by 
using the LRT formula Eq.(|l5|). 

The second step in the analysis is to determine the 
adiabatic V regime. This is done by checking whether 
the LRT calculation of 6E{t) is sensitive to cr. In order 
to get in the adiabatic regime LRT-based quantum cor- 
rections to the classical result, we have to take the level 
spacing statistics into account [|8[ We have pointed 
out the difficulty in observing such corrections. Rather 
we can improve over the LRT prediction by taking into 
account either higher order or Landau-Zener corrections 
to perturbation theory. 

The third step in the analysis is to calculate the para- 
metric perturbative profile P^'^{r). This is associated 
with getting a rough estimate for the core width F. Also 
the parametric scales Sxc and (5xp,t can be determined on 
the basis of this analysis. The former is determined from 
inspecting F, while the latter is determined by comparing 
6E of the LRT calculation to the P^'^ based prediction. 

The fourth step in the analysis is to distinguish be- 
tween the perturbative and the non-perturbative regime. 
For this purpose we have to look for the timing of the de- 
parture from parametric behavior, as in the lower panel 
of Fig. 4. For SE one can use the LRT based calculation 
of Eq.®. 

The fifth step in the analysis is to identify the sud- 
den regime. For this purpose we have to eliminate the 
" sudden time" from Eq. (^4|) . The problem is to estimate 
T{t) in the non-perturbative regime. If we have only the 
bandprofile as an input, then we can use the rough esti- 
mate T{t) ^ 6E{t). This is based on the assumption that 
in the non-perturbative regime the energy distribution is 
characterized by a single energy scale. 

The sixth step in the analysis is to make a sim- 
ulation of the parametric evolution. This is a rela- 
tively heavy task, but it is still much easier than mak- 
ing finite V simulation. [It requires merely diagonaliza- 
tions, while a temporal simulation requires an iterative 
procedure]. The main non-trivial effects that we have 
found in the analysis of the parametric evolution were: 



• Higher order tails grow up. 

• A non-perturbative core region develops. 

• The core-to-tail transitions are suppressed. 

• The tails are "pushed out". 

• Tail components are swallowed by the expanding 
core. 

The last two items are in the spirit of the "core-tail" 
theory, but go beyond the perturbative approximation 
Eq.(|lO|). Knowledge of the parametric evolution allows 
accurate determination of F, leading to a refined deter- 
mination of the V regimes. 

The seventh step in the analysis is to make simula- 
tion of the actual evolution. We would like to emphasize 
that this is the first time that the different dynamical 
scenarios (corresponding to the different V regimes) are 
illustrated in a numerical simulation. The only other nu- 
merical studies (mainly |l^) were too restricted in 
scope, and did not contain analysis of the stages in the 
evolution of the energy distribution. 

The only non-perturbative effect on the response that 
we have discussed so far is the Landau-Zener correction 
to the spreading in the adiabatic regime. Are there any 
other non-perturbative effects that affect the response? 
The immediate tendency is to regard LRT as the out- 
come of standard first order perturbation theory (Ap- 
pendi x D) . Then the question that arises is what happens 
if Eq. ( pDSp does not apply? 

Let us recall the answer in case of the parametric evo- 
lution of P{n\m). As 6x becomes larger we should be 
worried regarding the implications of having the effects 
that are listed at the end of Section 9. Having expand- 
ing "core" that "pushes out" the tails, and having grow- 
ing higher-order tail components, may sug gest that the 
first order calculation of the variance (Eq.(D7)) should 
be "corrected", and should include "higher order" terms. 
Does it mean that Eq. (|l^) underestimate the spreading? 
Or maybe we should assume that Eq.(|l^) provides the 
correct " trend" of higher order corrections? The tails in 
the vicinity of the core are growing slower, which can be 
regarded as a suppression of core-to-tail transitions due 
to the mixing. Consequently we would conclude that 
Eq.(|l6|) is an overestimation of spreading. Moreover, if 
we take Eq. (|l^) too seriously, beyond its regime of valid- 
ity, we would conclude that the spreading has saturation 
for large Sx. 

For the parametric evolution these possible specula- 
tions turn out to be wrong. The above effects are exactly 
balanced, and the LRT formula Eq.(p^ remains exact 
beyond any order of perturbation theory, which means 
that it is exact even in the non-perturbative regime where 
perturbation theory is not applicable. This claim has a 
simple derivation Note also that if the system has 
a classical limit, then the validity of Eq.(|l6|) can be es- 
tablished deep in the non-perturbative regime, where the 
semiclassical approximation becomes reliable. 

The question is whether this delicate balance is vio- 
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lated in case of finite V. For example, it may be that 
due to incomplete core-tail recurrences we shall have en- 
hanced spreading (compared with LRT). Unlike the para- 
metric case, we do not have a theoretical proof that ex- 
cludes such a possibility. In fact the contrary is the case. 
We have demonstrated that for an artificial (random ma- 
trix theory) model, the LRT formula cannot be trusted 
in the non-perturbative regime. Whether such effect is 
possible also for " quantized" system that possess a good 
classical limit has been left an open question. 

The possibility of having deviation from LRT in the 
non-perturbative regime was one important motivation 
for the present research. Clearly we did not witness such 
effect in our simulations (Fig. 8). This reflects that there 
is a clash between the semiclassical limit and the RMT 
limit. 



APPENDIX A: 
THE SPECTRAL FUNCTION Ft{uj) 

Given a function f{t') that describes the shape of the 
driving pulse during the time interval < t' < t, we 
define a spectral function Ft(w) by Eq.(^). This spec- 
tral function describes the spectral content of the driving 
pulse. Below we list some useful driving schemes. We use 
the notation 0(), which is defined by 8 (False) = and 
9 (True) = 1. 

For step function we have 

fit') = 9(0 < t') (Al) 
Fticu) = 1 (A2) 

For rectangular pulse we have 

fit') = 9(0 <t' <t) (A3) 
Fticj) - |l-e^-f (A4) 

If it is followed by a negative pulse pulse we get 

fit') = 9(0 < t' < t) - 9(| <t' <t) (A5) 
Ftiuj) = |1 - 26'^^* -(-e'^^'p iort>^ (A6) 



For linear driving we have 
fit') = t' 



Ftiuj) = i2(sinc(it^t))2 



(A7) 
(A8) 



where sinc(-) = sin(- )/(•). Note that for very large t we 
have 



Ftiuj) - t X 27r(5(w) 
Finally, for triangular pulse 

fit') = 2^9 (o < t' < I 



(A9) 



2(1-^)0(1 <t'<T) (AlO) 



we get 

Ftiuj) 

Ftiu) 



1 2 



T 

i^(sinc(icjt))^) for 0<i<- 



1 - 2c" 



for ^<t<T 



APPENDIX B: 
THE SPECTRAL FUNCTION Cico) 

Let us denote by F^n the matrix representation of 
some quantized observable FiQ,P), in the basis which 
is determined by some quantized chaotic Hamiltonian 
HiQjP). The bandprofile of Fm„ is conveniently char- 
acterized by the spectral function 



2n6 i - 1 (Bl) 



with implicit average over the reference state m. This 
is the power spectrum of the fluctuating quantity J^it), 
whose classical definition is J^it) = FiQit), Pit)), with a 
corresponding quantum-mechanical definition within the 
Heisenberg picture. The power spectrum of a fiuctuating 
quantity J- it) is defined as the Fourier transform of the 
corresponding correlation function C(r). 

Chaotic systems are characterized by fast decay of dy- 
namical correlations. We assume a separation of time 
scales between the (short) classical correlation time, and 
the (long) quantum-mechanical Heisenberg time. Thus, 
for times during which we can ignore the recurrences, we 
expect the following quantal-classical correspondence: 



Cir) 



C^\t) 



(B2) 



This correspondence implies that the envelope of Ciuj) is 
given by C'''(aj), and the following semiclassical expres- 
sion for the matrix elements follows [E4l: 



2TTh V ^ 



(B3) 



Taking into account the level spacing statistics, we de- 
duce the following relation p8[: 



Ciu^) 



Riuj) C'^iuj) 



(B4) 



where i?(w) oc {J2n '^('^ ~ ^nm))m is the two-point corre- 
lation function of the energy spectrum (Fourier transform 
of the spectral form factor). For the Gaussian orthogo- 
nal ensemble (GOE) of random matrix theory (RMT) 
the following result is well known: 



Rioj) = 1 - (sinc(7r?iw/A))^ 



(B5) 



and more generally it is common to assume Rioj) ^ ui^ 
for small ui. 
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APPENDIX C: CLASSICAL DERIVATION OF 
THE LRT FORMULA 

Consider ^{{Q, P;x{t)), and define the following time 
dependent quantities: 

m = -^{Q{t),p{t)-x{t)) (ci) 

E{t) = n{Q{t),P{t)-x{t)) (C2) 
E'{t) = n{Q{t),P{t)-xm (C3) 

Note that S{t) is the energy in the conventional sense, 
while £'{t) is the energy using a "fixed basis". We have 
the following relations 

d£{t) OH . . X X .X 



dS'jt) 
dt 



5x{t)T{t) 



(C5) 



The latter approximated equality strictly holds if the per- 
turbation 7i — Tio is linear with respect to the perturba- 
tion parameter 5x = x — Xq. From the equations above 
it follows that 



m - m 



x{t')T{t')dt' (C6) 



5x{t')T{t')dt' 



(C7) 



Obviously the two latter exp ress ions coincide for a cycle 
{x{t) — a:(0)). Squaring Eq.(|C6|), and performing a mi- 
crocanonical average over the (implicit) initial conditions 
(g(0),P(0)) one obtains 



5E{tf 



"'0 



f{t')f{t")C{t' - t")dt'dt" (C8) 



This can be written as Eq.(15). 



APPENDIX D: FIRST ORDER PERTURBATION 
THEORY AND LRT 

For convenience we assume that the perturbation is 
linear 'm 5x = x — xq. In complete analogy with the 
classical analysis we can work either using " fixed basis" , 
or else we can use the proper cc-dependent basis (the so 
called adiabatic representation). The respective matrix 
representations of the Hamiltonian are 



n ^^ + 5x{t)¥ 
H E + i(<)W 



(Dl) 
(D2) 



where E is a diagonal matrix. Note that in the first 
equation E and F are calculated for x — xq, while in 
the second equation there is an implicit x{t) dependence. 
The matrix elements of F arc 



m 

dx 



(D3) 



The off diagonal matrix elements of W are 



nm — - ( ^ 



dt 



En—E„ 



(D4) 



and we use the 'gauge' convention W„m = for n=m. 
(Only one parameter is being changed and therefore 
Berry's phase is not an issue). The derivation of Eq.(D2) 
is standard, and can be found in Section 11 of [l^. 



Using first order perturbation theory with Eq.(D2) we 



get 



Pt{n\m)^5nm + A^Ft 



w 



(D5) 



Note that for a cycle {f{t) = /(O) = 0), the same re- 
sult is o btained via first order perturba tion theory with 
Eq. ( D 1 ) . As a global approximation Eq. ( D5 ) is valid only 
in the standard perturbative regime. In the extended 
perturbative regime it is valid only for the first order tail 
region (see Section 8). The expression for the first order 
tail region can be written in a concise way as 



Ann 



1 



flUJ„ 



(D6) 



The corresponding global approximation is given by 
Eq.(|0l). 

Assuming that the variance is dominated by the first 
order tail component, we get the LRT result 

SE{t)^ = J2 Pt{n\m) X (^„ - E^)^ 



A' X 



^F,{u;)C{u;) (D7) 



An implicit average over m is assumed. Note that the 
latter formula does not contain h. This "restricted" 
quantal-classical correspondence holds only for the vari- 
ance. Higher moments of the energy distribution are typ- 
ically much smaller compared with the classical expecta- 
tion, and scale like h to the power of the moment order 
minus 2. 
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Fig.l. Images of time evolution within one-pulse period for V = 100 
(upper panel) and for V = 1 (lower panel). Each column is a profile 
of the probability distribution Pt{r) for a different time step t. The 
first 7 time steps are log spaced, while the rest are linearly spaced. 
The V = 100 evolution is predominantly parametric. 
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Fig. 2. Spreading profile at the end of one pulse period for V = 100 
(upper panel) and for V = 10 (lower panel). The thin solid lines are 



Eq.(0) 



with r = 0. The dashed lines are for Gaussian with the 



same variance. 
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Fig. 3. The width r{t) as a function of time. The numerical de- 
termination of r is explained after Eq. (^. The horizontal axis is 
the scaled time Vt. The different curves correspond to the differ- 
ent velocities V = 100, 80, 40, 20, 10, 7, 5, 3, 1, 0.5, 0.2, 0.1, 0.07, 0.05. 
The solid lines highlight the velocities V = 100 (upper most), and 
V = 10, and V = I, and V = 0.1. The lower panel is the same 
as in the upper figure, but in log-log scale, and only half period is 
displayed. 
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Fig. 4. The spreading 5E(t) of Eq.(|g) as a function of time, for 
different V simulations. The upper panel is in one to one cor- 
respondence with the upper panel of the previous figure. In the 
lower panel the data of the upper panel are presented in a differ- 
ent way: the vertical axis is the relative spreading; The horizontal 
axis is log scale; And only half period is displayed. The location 
of (Sxprt = 0.05 is indicated by vertical line. See the text for more 
details. 
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Fig. 5. The band profile, as defined in Appendix B. The two curves 
are the outcome of two different numerical procedures: The thin 
line is based on direct evaluation of matrix elements, while the thick 
line is deduced form the evolution over an infinitesimal time step 
(see Section 12 for further details). A third possi ble procedure (not 
displayed) is to use the semiclassical recipe Eg. ( [B3| ) 
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Fig. 7. LRT calculation of the spreading versus time, using Eq.(p^ 
with the bandprofile as an input. The horizontal axis is Vt. This 
Figure is in one to one correspondence with the upper panel of 
Fig. 4. One observes that for small V the calculation underestimates 
the observed spreading. See discussion in Section 12. 




10"^ 10"^ 10"' 

5x 



Fig. 6. Parametric evolution: The horizontal axis is 5x either in 
linear scale (upper panel) or in logarithmic scale (lower panel). 
The dotted line is the width F from (M). The thick line is 5E from 
(^), while the dashed line is the LRTestimate Fg l]lp[ ). The thin 
line is the perturbative estimate Eg. (H) using Eg. (llHj! From this 
plots we can determine 5xc = 0.006 and Sxpit = 0.05. 
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Fig. 8. upper panel: Response for one-pulse versus V. The Thick 
line is 5E at the end of the pulse. The thiiL(solid and dotted) lines 
correspond to the LRT calculation Eg. (ha). For the calculation 
of the lower thin lines we have set the near-level coupling a = 0. 
The dashed line is the spreading at the end of half pulse period. 
The dash-dot line is the slope that corresponds to Landau-Zmw 
spreading mechanism. The lower panel is the difference Eq.(|l3|) 
from Gaussian line shape. 



